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After a short excursion from discovery of Brownian motion to the Richardson " law 
of four thirds" in turbulent diffusion, the article introduces the Levy flight superdif- 
fusion as a self-similar Levy process. The condition of self-similarity converts the 
infinitely divisible characteristic function of the Levy process into a stable char- 
acteristic function of the Levy motion. The Levy motion generalizes the Brown- 
ian motion on the base of the a-stable distributions theory and fractional order 
derivatives. The further development of the idea lies on the generalization of the 
Langevin equation with a non-Gaussian white noise source and the use of func- 
tional approach. This leads to the Kolmogorov's equation for arbitrary Markovian 
processes. As particular case we obtain the fractional Fokkcr-Planck equation for 
Levy flights. Some results concerning stationary probability distributions of Levy 
motion in symmetric smooth monostable potentials, and a general expression to 
calculate the nonlinear relaxation time in barrier crossing problems are derived. 
Finally we discuss results on the same characteristics and barrier crossing problems 
with Levy flights, recently obtained with different approaches. 
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ties of a medium, strengthened with increas- 
ing medium temperature, with a diminution 
of its viscosity and sizes of particles. But R. 
Brown considered the motion of the parti- 
cles (not being atoms of course) as their own 
property and said nothing about atoms or 
molecules. 

It should pass almost 8 decades before two 
physicists Albert Einstein [Einstein, 1905] 
and Marian von Smoluchowski [Smolu- 
chovski, 1906] found the physical explanation 
of Brownian motion. It was based on consid- 
eration of thermal motion of molecules sur- 
rounding the Brownian particle. The history 
of the further study of Brownian motion is as- 
sociated with names of Langevin [Langevin, 
1908], Perrin [Perrin, 1908], Fokker [Fokker, 
1914], Planck [Planck, 1917], Uhlenbeck, 
Ornstein [Uhlenbeck & Ornstein, 1930], 
Chandrasekhar [Chandrasekhar, 1943] and 
other well-known physicists. However, for 
the first time the diffusion equation appeared 
in the thesis of Louis Bachelier [Bachelier, 
1900], a student of A. Poincare. His thesis, 
entitled " The theory of speculations" , was de- 
voted to the study of random processes in 
market prices evolution. 

It is astonishing, how the same diffusion 
equation can describe the behavior of neu- 
trons in a nuclear reactor, the light in at- 
mosphere, the stock market values rate on fi- 
nancial exchange, particles of flower dust sus- 
pended in a fluid and so on. The fact that 
completely different by nature phenomena 
are described by identical equations is a di- 
rect indication that the matter concerns not 
the concrete mechanism of the phenomenon, 
but rather the same common quality of whole 
class of similar phenomena. The statement 
of this quality in terms of physical laws and 
mathematical postulates or definitions allows 
to liberate a given pattern from details, which 
are not influencing essentially the physical 
process, and to explore the obtained model 
through general laws. This is a typical situa- 
tion for statistical physics and applied math- 
ematics. The new approaches proposed by 



Einstein, Smoluchowski and Langevin to de- 
scribe the Brownian motion, in fact, open 
the door to model a great variety of natu- 
ral phenomena. At the same time for mathe- 
maticians, whose achievements built the the- 
ory of random processes, the first object of 
its application became the Brownian motion. 
The major contribution to the mathemat- 
ical theory of Brownian motion has been 
brought by N. Wiener [Wiener, 1930], who 
has proved that the trajectories of Brown- 
ian process almost everywhere are continuous 
but are not differentiable anywhere. Along 
with Wiener the mathematical aspects of 
Brownian motion were treated by Markov, 
Doob, Kag, Feller, Bernstein, Levy, Kol- 
mogorov, Stratonovich, Ito and others [Doob, 
1953; Kac, 1957; Feller, 1971; Levy, 1925, 
1965; Kolmogorov, 1941; Stratonovich, 1963, 
1967, 1992; Ito, 1944, 1946, 1965]. 

Two important properties are intrinsic to 
the homogeneous Brownian motion: the dif- 
fusion packet initially concentrated at a point 
takes later the Gaussian form, whose width 
grows in time as t 1 / 2 . This kind of diffusion 
was called the normal diffusion. 

Twenty years later Einstein, Smolu- 
chowski and Langevin works, L. Richard- 
son published the article [Richardson, 1926] 
where he presented empirical data being in 
contradiction with the normal diffusion: the 
size A of an admixture cloud in a turbu- 
lent atmosphere grows in time proportionally 
to i 3 / 2 , that is much faster then in the nor- 
mal case (t 1 ' 2 ). This turbulent diffusion was 
the first example of super diffusion processes, 
when A oc t 7 with 7 > 1/2. 

The phenomenon has been interpreted 
as a diffusion process with a variable dif- 
fusivity D{r) oc r 4 / 3 . This Richardson's 
"law of four thirds" was grounded theo- 
retically by Russian mathematicians A. N. 
Kolmogorov [Kolmogorov, 1941] and A. M. 
Obukhov [Obukhov, 1941] as a consequence 
of the self-similarity hypothesis of locally 
isotropic turbulence. However, the fact that 
the diffusivity should depend not on the co- 
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ordinates (turbulent medium is supposed to 
be homogeneous in average), but on the scale 
or distance I between a pair of diffusing parti- 
cles, creates essential difficulties both to find 
a solution to the equation and for its inter- 
pretation. 

As Monin showed in Ref. [Monin, 1955], 
the same law of the diffusive packet widen- 
ing with time can be obtained in the frame- 
work of the homogeneous Markovian pro- 
cesses family, when the characteristic func- 
tion of the spatial distribution of the diffusive 
substance 

P(k,t) = U ikX{t) ) = / e ikx P{x,t)dx 

* J -co 

(1) 

obeys the equation 

= -c\k\°p{k,t) (2) 

with C being a positive constant. Under ini- 
tial condition P(x, 0) = 5(x) we obtain from 
Eq. © 

P(k,t) = exp{-C\k\ a t}, (3) 

where a = 2/3. This is nothing but the 
characteristic function of the a-stable Levy 
distribution, and the random process itself 
is the Levy motion (Levy flights). Later, 
anomalous diffusion in the form of Levy 
flights has been discovered in many other 
physical, chemical, biological, and finan- 
cial systems [Shlesinger et al, 1993; Met- 
zler & Klafter, 2000, 2004; Metzler et al., 
2007; Brockmann & Sokolov, 2002; Eliazar & 
Klafter, 2003; Barkai, 2004; Chechkin et al, 
2006; Mandelbrot, 1997; Mantegna, 1991]. 

Levy flights are stochastic processes char- 
acterized by the occurrence of extremely long 
jumps, so that their trajectories are not con- 
tinuous anymore. The length of these jumps 
is distributed according to a Levy stable 
statistics with a power law tail and diver- 
gence of the second moment. This pecu- 
liar property strongly contradicts the ordi- 
nary Brownian motion, for which all the mo- 
ments of the particle coordinate are finite. 



The presence of anomalous diffusion can be 
explained deviation of the real statistics 
of fluctuations from the Gaussian law, giv- 
ing rise to the generalization of the central 
limit theorem by Levy and Gnedenko [Levy, 
1925, 1937; Gnedenko & Kolmogorov, 1954]. 
The divergence of the Levy flights variance 
poses some problems as regards to the physi- 
cal meaning of these processes. However, re- 
cently the relevance of the Levy motion ap- 
peared in many physical, natural and social 
complex systems. The Levy type statistics, 
in fact, is observed in various scientific areas, 
where scale invariance phenomena take place 
or can be suspected. Among many interest- 
ing examples we cite here chaotic dynamics of 
complex systems [Zaslavsky, 2005; Solomon 
et al., 1993, 1994], diffusion and annihilation 
reactions of Levy flights with bounded long- 
range hoppings [Albano, 1991], front dynam- 
ics in reaction-diffusion systems with Levy 
flights [del-Castillo-Negrete et al, 2003], frac- 
tional diffusion [West et al, 1997; Chaves, 
1998], thermodynamics of anomalous diffu- 
sion [Zanette & Alemany, 1995], dynamical 
foundation on noncanonical equilibrium [An- 
nunziato et al., 2001], quantum fractional ki- 
netics [Kusnezov et al, 1999], trapping dif- 
fusion [Vazquez et al, 1999] , Levy diffusion 
processes as a macroscopic manifestation of 
randomness [Grigolini et al, 1999; Bologna 
et al, 1999], diffusion by flows in porous 
media [Painter, 1996], two-dimensional Levy 
flights [Desbois, 1992], persistent Levy mo- 
tion [Chechkin & Gonchar, 2000], self- 
avoiding Levy flights [Grassberger, 1985], 
Levy flights with quenched noise ampli- 
tudes [Kutner & Maass, 1998], cooling down 
Levy flights [Pavlyukevich, 2007], branching 
annihilating Levy flights [Albano, 1996], ran- 
dom Levy flights in the kinetic Ising and 
spherical models [Bergersen & Racz, 1991; 
Xu et al, 1993], plane rotator in presence of 
a Levy random torque [Caceres, 1999], fluc- 
tuations and transport in plasmas [Chechkin 
et al., 2002b; Lynch et al., 2003], transport in 
stochastic magnetic fields [Zimbardo & Vel- 
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tri, 1995], Levy flights in the Landau- Teller 
model of molecular collisions [Carati et al, 
2003], subrecoil laser cooling [Bardou et al, 
1994, 2002; Reichel et al, 1995; Schaufler et 
al, 1997, 1999], scintillations and Levy flights 
through the interstellar medium [Boldyrev 
& Gwinn, 2003], Levy flights in cosmic 
rays [Wilk & Wlodarczyk, 1999], anoma- 
lous diffusion in the stratosphere [Seo & 
Bowman, 2000], long paleoclimatic time se- 
ries of the Greenland ice core measure- 
ments [Ditlevsen, 1999a], seismic series and 
earthquakes [Posadas et al, 2002; Sotolongo- 
Costa et al, 2000], signal processing [Nikias 
& Shao, 1995], time series statistical analysis 
of DNA [Scafetta et al, 2002], primary se- 
quences of proteinlike copolymers [Govorun 
et al, 2001], spatial gazing patterns of bac- 
teria [Levandowsky et al, 1997], Levy-flight 
spreading of epidemic processes [Janssen et 
al, 1999], contaminant migration by biotur- 
bation [Reible & Mohanty, 2002], flights of 
an albatross [Viswanathan et al, 1996; Ed- 
wards et al, 2007], fractal time in animal 
behaviour of Drosophila and animal locomo- 
tion [Cole, 1995, Seuront et al, 2007], finan- 
cial time series [Mandelbrot, 1963; Bouchaud 
& Sornette, 1994; Mantegna & Stanley, 1996, 
1998; Chowdhury & Stauffer, 1999], human 
stick balancing and Levy flights [Cabrera & 
Milton, 2004] and human memory retrieval 
as Levy foraging [Rhodes & Turvey, 2007]. 
Experimental evidence of Levy processes was 
also observed in the motion of single ion 
in one-dimensional optical lattice [Katori et 
al, 1997] and in the particle evolution along 
polymer chains [Sokolov et al, 1997; Lomholt 
et al, 2005], and in self-diffusion in systems 
of polymer like breakable micelles [Ott et al, 
1990]. 

Levy flights special class of Marko- 

vian processes, therefore the powerful meth- 
ods of the Markovian analysis are in force in 
this case. We mean a possibility to investi- 
gate the stationary probability distributions 
of superdiffusion, the first passage time and 
the residence time characteristics, the spec- 



tral characteristics of stationary motion, etc. 
Of course, this type of diffusion has a lot 
of peculiarities different from those observed 
in normal Brownian motion. The main dif- 
ference from ordinary diffusion consists in 
replacing the white Gaussian noise source 
in the underlying Langevin equation with a 
Levy stable noise. 

In the first part of the present paper, we 
give a short introduction to the Levy motion. 
Being a generalization of the Brownian diffu- 
sion, it takes an intermediate place between 
Brownian motion and Levy processes (i.e. 
infinitely divisible processes, see [Bertoin, 
1996]) in the random process hierarchy sys- 
tem. The Levy motion is introduced as a 
self-similar Levy process. 

The second part of this paper is devoted 
to the stationary probabilistic characteristics 
and the problem of barrier crossing for Levy 
flights. We use functional approach to derive 
the generalized Kolmogorov equation directly 
from Langevin equation with the Levy pro- 
cess [Dubkov & Spagnolo, 2005]. In particu- 
lar case of Levy stable noise source we obtain 
the Fokker-Planck equation with fractional 
space derivative. Starting from this equation 
we find the exact stationary probability dis- 
tribution (SPD) of fast diffusion in symmet- 
ric smooth monostable potentials for the case 
of Cauchy stable noise. Specifically, we con- 
sider potential profiles U(x) = / yx 2m /(2m), 
with odd and even m, useful to describe the 
dynamics of overdamped anharmonic oscil- 
lator driven by Levy noise. We find that 
for Levy flights in steep potential well, with 
potential exponent 2m greater or equal to 
four, the variance of the particle coordinate 
is finite [Dubkov & Spagnolo, 2007]. This 
gives rise to a confined superdiffused mo- 
tion, characterized by a bimodal stationary 
probability density, as previously reported in 
Refs. [Chechkin et al, 2002a, 2003a, 2004, 
2006]. Here we analyze the SPD as a func- 
tion of a dimensionless parameter f3, which is 
the ratio between the noise intensity D and 
the steepness 7 of the potential profile. We 
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find that the SPDs remain bimodal with in- 
creasing j3 parameter, that is with decreasing 
the steepness 7 of the potential profile, or by 
increasing the noise intensity D. 

The particle escape from a metastable 
state, and the first passage time density 
have been recently analyzed for Levy flights 
in Refs. [Ditlevsen, 1999b; Rangarajan & 
Ding, 2000a, 2000b; Buldyrev et ai, 2001; 
Chechkin et ai, 2003b, 2005, 2006, 2007; Dy- 
biec & Gudowska-Nowak, 2004; Dybiec et ai, 
2006, 2007; Bao et ai, 2005; Ferraro & Za- 
ninetti, 2006; Imkeller & Pavlyukevich, 2006; 
Imkeller et ai, 2007; Koren et ai, 2007]. The 
main focus in these papers is to understand 
how the barrier crossing behavior, according 
to the Kramers law [Kramers, 1940], is mod- 
ified by the presence of the Levy noise. Fi- 
nally we discuss briefly some results on the 
barrier crossing events with Levy flights, re- 
cently obtained with different approaches. 

II. LEVY PROCESSES 

To see better a place of the diffusion pro- 
cesses under consideration among other ran- 
dom processes we shall remind some defi- 
nitions. We restrict ourselves to the one- 
dimensional case when X, x G (—00,00) and 
t > 0. 

A random process {X(t),t > 0} is a set of 
random variables X (t), given on the same 
probability space and corresponding to any 
possible time t > 0. 

A random process {X(t),t > 0} is called 
a Markovian process, if for any n > 1 and 
h <t 2 < ... <t n <t P(X(t) < x \Xfa) = 

X 1 ,...,X(t n ) = X n ) = P(X(t) < x\X(t n ) = 

x n ). The Markovian property is interpreted 
as independence of future from the past for 
the known present. P. Levy states this prop- 
erty by the sentence "the past influences the 
future only through the present" and under- 
lines analogy to Huygens' principle ("it is 
possible to say, that it is Huygens ' principle 
in calculation of probabilities" [Levy, 1965]). 
A random process {X(t),t E T} is called 



the process with independent increments if for 
any n > 1 and t\ < ti < . . . < t n < t random 
variables X (0) , X fa) - X (0) , . . . , X (t n ) - 
X (t„_i) are mutually independent. The ran- 
dom variable X (0) is called the initial state 
(value) of the process, and its probability dis- 
tribution is called initial distribution of the 
process. P. Levy named such processes addi- 
tive. Obviously, they belong to the class of 
Markovian processes. 

A random process with independent incre- 
ments is called homogeneous or stationary, if 
the random variables X (t + r) — X (t) have 
distributions which are independent on t: 

P{X(t + r)-X(t) <x} = F(x,r). (4) 

P. Levy named such processes linear, re- 
marking that among them "there are also 
processes distinct from Brownian" . Now, 
Bertoin [Bertoin, 1996] and Sato [Sato, 1999] 
use the term Levy processes for the processes 
with stationary independent increments. 

One can paraphrase the definition by say- 
ing that {X(t),t > 0} is a Levy process if, for 
every t, r > 0, the increment X(t + r) —X(t) 
is independent on the process {X(t'), < t' < 
t} and has the same law as X(t). In partic- 
ular, X(0) = 0. 

We will denote the Levy process L(t). As 
it follows from the evident decomposition 

the random variable L(t) can be divided into 
the sum of an arbitrary number of indepen- 
dent and identically distributed random vari- 
ables. In other words, the probability distri- 
bution of L(t) belongs to the class of infinitely 
divisible distributions [de Finetti, 1929, 1975; 
Khintchine, 1938; Khintchine & Levy, 1936; 
Levy, 1937; Gnedenko & Kolmogorov, 1954; 
Feller, 1971; Mainardi & Rogosin, 2006]. 
Hence, we can express the second characteris- 
tics, i.e. the logarithm of characteristic func- 
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tion of the random variable L(t) in the Levy- 
Khinchine form [Feller, 1971] 



0(M) =lnP{k,t) = ln (e lkL{t) ) 



+00 



(e ikx 



1 ■ \ POM) , ta\ 
tksmx) — ax, (b) 



x- 



where p(x, t) is the canonical measure density 
(with respect to the first argument). 

For two consecutive non-overlapping time 
intervals t x and t 2 we have 



L{t x + t 2 )±L(t x )+Uf2) 



(7) 



where L(ti) and L(t 2 ) are mutually inde- 
pendent random variables and the symbol = 
means the equality of distributions of the cor- 
responding random variables. Therefore, 

P(k,t 1 +t 2 ) = P(k,t 1 )P(k,t 2 ) (8) 

or 

0(Mi + *2) = 0(Mi) + #(M2). (9) 

According to Eqs. ([6]) and (Q we have 

p(x,tx + t 2 ) = p{x,ti) + p(x,t 2 ). (10) 

The differentiable solution of Eq. (TlOl . re- 
garding t, is only linear one 

p{x,t)=tp{x). (11) 

So, from Eq. §b§ we obtain 



•{k,t)=t f °° (e ikx -1-ik 



sin 2) — 75-dx, 
(12) 

where the kernel p[x) > 0. Note that the 
last term in the bracket, — ik sin x, serves to 
ensure the convergence of the integral and 
can be omitted if the integral converges it- 
self. Choosing 

p(x)=6(x), (13) 
and taking into account that 



e — 1 — ik sinx = —k x /2 + o(x 



0. 



we arrive at the normalized Brownian mo- 
tion B(t) (Wiener process) with characteris- 
tic function 

P(k,t) =exp{-tk 2 /2}. (15) 



III. SELF-SIMILARITY (SCALING) 

The self-similarity (scaling is a synonym 
of self- similarity) is a special symmetry of a 
system (process) revealing that the modifica- 
tion of the scales of one variable can be com- 
pensated by the homothetic transformation 
of the others. For example, if the state of a 
system is characterized by function u(x,t), 
where x is the coordinate, t is the time, 
the requirement of invariance with respect 
to scale transformations x — > x' = kx and 
t — > t' = It, looks like 



u(x,t) = k a l 5 u (kx,lt) 



(16) 



where k and I are positive, and a and S are ar- 
bitrary numbers. By choosing k a = I = m/t, 
where m > is a parameter of similarity, 
we obtain a self-similar form for the function 

u(x, t) 

u (x, t) = (m/t) 1+S u ( K [m/t) 1 ^ a x, m \ . (17) 

In our case such a function is the probability 
density function P(x,t). The normalization 
condition 



/+00 
P{x,t) dx = 1 
-00 



(18) 



and the principle of self-similarity (fTTj) give 
1 + 5 = 1/a and lead to the representation 
(m = 1) 

P(x,t) =r 1/a g {a) (xr 1/a ) , (19) 



where 

gW{x)=P(x,l). (20) 

It must be emphasized that we can define the 
self-similarity of a random process X(t) with 
stationary increments as usually (see, for ex- 
ample, Ref. [Chechkin et al., 2002c]) 

(14) X (t + kt) - X it) = K H [X(t + r)-X (t)] . 



In such H = 1/a. In terms of charac- 

teristic functions we have 



P(k,t) =g {a) (kt l/a ) 



(21) 
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with g( a > obeying the equation 



(k(t 1 + t 2 f«)=gM (kt{ /a )~g^ {ktf 

(22) 

which follows from Eq. ([S]). 

Let be a random variable described 
by the characteristic function 

(k) = (exp {ikY {a) }) . (23) 

Obviously, 



~(a) f kt V*\ = ^ exp {ikY (a) t 1/a }) (24) 

determines the random variable t^^Y^ sat- 
isfying the relation 

( tl + t2 f« y(«)4 t ;/ a yW + tl /a Y 2 {a \ (25) 

where and Kj are independent copies 
of random variable Y^. This relation is a 
definition property of strictly stable random 
variables with a characteristic index a. We 
arrive, therefore, at the very important sub- 
family of the Levy processes called the Levy 
motion (often, the term " Levy flights" is used 
clS cL synonym) . 

IV. STABLE RANDOM VARIABLES 

To find an explicit expression for the sta- 
ble characteristic functions, we can proceed 
by two equivalent ways: (i) using the general 
representation of infinitely divisible charac- 
teristic functions (|T2|) or (ii) using the stabil- 
ity property (12"2"I) . We choose the latter way. 

Let us introduce the second characteristic 



</> (a) (k) =\ng {a) (jfc) 



(26) 



so the property ff22l) of a strict stability takes 
the form 

^ (Q) (AiJfc) + tp {a) (\ 2 k) = ^ (Q) (AJfc) , (27) 
where 



Extending this relation to the sum of ar- 
bitrary number n of identically distributed 
(Ai = A 2 = • • • = A n = 1) terms, we obtain 

(Jfe) = ^ (n 1/a k) . (29) 

According to the property 



^(°o (-k) = [v> (Q ° (fc) 



(30) 



it is enough to determine the function ip^ (k) 
for positive arguments k > 0. Taking into ac- 
count its continuity in a neighborhood of the 
origin and the initial condition of the charac- 
teristic function 



i) {a) (0) = 0, 



(31) 



we obtain that 



(jfe) = const ■ k a (k > 0, a > 0) 



and 



¥ a) (k) = -k a (cb-id), 



(32) 
(33) 



where Co and C\ are arbitrary real constants. 
Since the characteristic function satisfies the 
requirement 



~g^\k) 



< 1, 



then 



Re{^ {a) {k)} < 



(34) 



(35) 



and the real constant Co should be positive. 
On the other hand, from Eqs. (|23|) and (|26|) 
we have 

£(«)'(o) = # (a) "(o) = -([r (Q) ] 2 ) 

(36) 

and 

^(«)"(o) = -([y (a) ] 2 ) + (r (Q) ) 2 

= -Var{F (Q) } < 0. (37) 

Calculating the second derivative from 
Eq. fl33|) . we obtain 



A = (A? + A£) 1/q . 



(28) V (Q),/ (^) = ~(co - ic x )a(a - l)*; " 2 . (38) 
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As one can see from Eq. fl38l) . we have: (i) for 
a = 2 the variance is finite (thus the constant 
ci should be equal to zero since the variance 
is real); (ii) for a < 2 and k — > we obtain 
the infinite variance (in this case c\ does not 
play any role), and (iii) for a > 2 and k — ► 
the derivative gives zero. This means that in 
the expression for the second moment 

/+oo 
x 2 g ia \x)dx, 
-00 

the function g( a )(x) ceases to be a probabil- 
ity density, when the characteristic index ex- 



ceeds the boundary value a = 2. We come 
to the conclusion that a range of values of 
the characteristic index a is the interval (0, 2] 
closed on the right. Because of Cq > and 
— oo < Ci < +oo, the constants of Eq. f[3"3"j) 
can be put in the form 

c = l, c 1 = /3tg (ott/2), -1 < /? < 1. 

(39) 

Therefore, the characteristic function 
g( a 'P>(k) of the strictly stable probability 
distribution g( a >P'(x), with parameters a and 
/5, is given by the formula 



~g^\k) 



exp 



/a;7r\ ' 
1 — ip tan I — 1 sgn k 



(40) 



where sgn x is the sign function. The charac- 
teristic index a (with a < 2) determines the 
decreasing rate of the large values probability 
for stable distributions 



p{|y(a,/?)| > A | x A -a A 



oo. 



(41) 



The parameter /3 characterizes the asymme- 
try of the distributions: for (3 = the sta- 
ble distribution is symmetric. The class of 
the symmetric stable distributions, besides 
the above-mentioned Gaussian distribution, 
includes also the Cauchy distribution 



9 {lfi) (*) 



7T(1 +X 2 ) 



with the characteristic function 



~g {1 ' 0) (*) 



exp{- |fc|}. 



(42) 



(43) 



For a < 1 the distributions with extreme 
values of asymmetry (3 are located on a semi- 
axes only: positive if (3 = 1 or negative if 
/3 = — 1. One of these well-known one- side 
distribution is the Levy - Smirnov distribu- 
tion 



x 3//2 exp 



1 

"2x 



x>0. 
(44) 



The detailed exposition of properties of 
stable random variables and their distribu- 
tions can be found in the books [Khintchine, 
1938; Levy, 1965; Gnedenko & Kolmogorov, 
1954; Feller, 1971; Bertoin, 1996; Sato, 1999; 
Uchaikin & Zolotarev, 1999b]. We shall un- 
derline here only the fact that all members 
of the set of stable distributions are charac- 
terized by presence of "heavy" (power- type) 
tails and, as a consequence, of infinite vari- 
ance, and that concerns all of them, except 
the Gaussian (normal) distribution. From 
the point of view of the whole " noble family" , 
the Gaussian distribution should be looking 
defiantly abnormal, monstrous, ugly duckling 
among white swans. For us (at least, for 
many of us) the infinite variance associates 
with an infinite error (what is not the truth), 
with an infinite energy or with something else 
what does obviously not have any physical 
sense. 

Professional physicists will recognize in 
Cauchy density a natural profile of a radia- 
tion line or the cross-section formula for reso- 
nances in nuclear reactions, and they will re- 
member the Holtzmark distribution describ- 
ing fluctuations of electric field strength cre- 
ated by Poisson ensemble of point-like ions 
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in plasma and the fluctuations of the gravita- 
tion field of stellar systems. But at this stage 
their acquaintance with the stable laws usu- 
ally comes to the end. This should be caused 
by the circumstance that the stable densi- 
ties, as a rule, are not expressed in terms of 



elementary functions: the above mentioned 
formulas fully exhaust a set of "convenient" 
distributions. There are some more distribu- 
tions which are expressed through the known 
special functions, like the following one (see 
Ref. [Garoni & Frankel, 2002]) 



^(2A0) = 

V ' 2^ ' 

where W^ u (x) is the hypergeometric Whit- 
taker function. However, in the age of com- 
puters the lack of simple formulas has no so 
strong importance. In fact, if simple or com- 
plex formulas can be processed by a com- 
puter, never mind, it computes fast and very 
well, provided that we check the full process. 

The properties of "anomalous" stable dis- 
tributions are really remarkable. If we shall 
sum up n independent random variables, 
distributed under the same stable law, the 
breadth of this new distribution will grow 
proportionally to n 1 ^, and the breadth of 
distribution of arithmetic mean will grow as 
n i/a-i p Qr a _ i w hen not only the vari- 
ance diverges, but also the expectation value 
does not exist, the width of arithmetic mean 
distribution remains constant! And if the re- 
sults of your measuring are distributed by 
the Cauchy law, the repetition of measuring 
cannot decrease the "statistical error" in no 
way. If a < 1, increasing the number of terms 
in the sum of results in widening of " sample 
mean" distribution! Clearly, the law of large 
numbers does not work here because the ex- 
pectation value does not exist. 



In summary we remark that the whole set 
of stable laws appear as limiting distributions 
in the generalized central limiting theorem: 
any other laws cannot be limiting ones. This 
is, certainly, their most important advantage. 



Tf 2 ) (2k 2 ) • (45) 

I 

V. STABLE PROCESSES AND LEVY 
MOTION 

Having designated the random realiza- 
tion of the processes under consideration 
as lS a, P\t), we write the condition of self- 
similarity as 

£(°./3)( f ) =t V<»yM ) (46) 

where Y^ a '^ is the strictly stable random 
variable, with parameters a and (5. The set 
of such processes is sometimes called stable 
processes. 

The random process {X (t) ,t e T} is 
called stable (strictly stable), if all its finite- 
dimensional distributions are stable (strictly 
stable). This definition generalizes the con- 
cept of Gaussian process, not restricted by 
the requirements of homogeneity and inde- 
pendence of increments. Thus, however, it is 
necessary to introduce the concept of mul- 
tivariate stable distribution or multivariate 
stable vector. 

The random vector Y = (Y" 1? . . . , Y m ) is 
called stable random vector in K m , if for 
any positive numbers Ai,A 2 there are posi- 
tive number A and vector c e SR m , such that 

AiY' + A 2 Y" = AY + c, (47) 

where Y',Y" are independent copies of the 
random vector Y. 

The stable random vector Y is called 
strictly stable if the last equality, with c = 0, 
holds true for any A x and A 2 . 
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The stable random vector Y is called sym- 
metric stable random vector, if it satisfies the 
relation 



P{Y e A} = P{-Y G A} 



(4£ 



for any Borel set A C 3ft m . Similarly to the 
one-dimensional case, the symmetric vector 
is strictly stable (the inverse statement, cer- 
tainly, is not true). 

Recall the standard definition. The ran- 
dom process {L^(t),t> 0} is called (stan- 
dard) a -stable Levy-motion with parameters 
< a < 2, -1 < < 1, if 

I ]Jpt,p) (g) = almost certainly; 

2. {l^^ (t) ,t > o| is a process with in- 
dependent increments; 



r>,/?) (t + r)- (t) 
at any t and r. 



. T l/ay(a,P) 



For the sake of brevity we shall call it L^ a '^'- 
process, then the Wiener process will be des- 
ignated as L( 2,0 )-process. 

For better understanding of the principal 
difference between L-processes with a = 2 
and a < 2, consider the third property for 
the Wiener process, namely, the Lindeberg 
condition which reflects the continuity of its 
trajectories. According to the continuity cri- 
terion of a random process [Loeve, 1963], for 
a = 2 and r — > we have 



p|| L (2,o) (t + r) _ L (2,o) W | > A | P j|y(2,o)| > a/v^F} i 



e- z2/4 dz. (49) 



Evaluating the indeterminate form by From Eq. (14ip and using the property 3, for 
l'Hopital's rule, we obtain for r —>■ a < 2 and r — > 0, we have 

-1= 4- r e~^dz = -^-r-^e- A2 ^ - 0. 

(50) 



P { |L^) (t + r) - (t) | > A} P {r 1 /" | Y | > A} 



p||y(«,/3)| > Ar" 1 /"} 



oc A~ a > 0. 



(51) 



Thus, L( 2, °) is the only L^ a '^ process pos- index a is the fractal dimensionality of the 
sessing continuous trajectories. As it was Levy process trajectories, 
shown in Ref. [Seshadri & West, 1982], Levy The size of a diffusion package for a < 2 
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grows with time faster than y/i, namely, pro- under initial condition 
portionally to t 1//o , and its shape differs from 

the Gaussian law. The variance is infinite, P^ 2,0 \k, 0) = 1. (54) 
but nothing interferes in using any other mea- 
sure of width, for example, the width on the Factor -k 2 is the Fourier image of the one- 
half of peak height, or the width of the in- dimensional Laplace operator Ai = d 2 /dx 2 . 
terval containing some fixed probability. For The inverse transformation yields the partial 
example, if we consider the fractal moment differential equation 
of Levy process increments 

dP™(x,t) _ &P™(x,t) 
L^(t) -L^)(0)| ) = / P(x,t)\x\ 5 dx, dt ~ dx~ 2 ' 1 j 

' with initial condition 

which is finite for < 5 < a, we immediately 

obtain from Eq. (19) p( ( x > °) = S ( x )- (56) 

jja,p) ^ _ jjufi) ^ 5 \ ^ ^s/a For the symmetric Levy motion with an arbi- 

/ trary a, the corresponding expression of the 
characteristic function reads 

and, as a result (see Ref. [Chechkin et al, 

2002c]) P<«.°>(M) = e-«l*l a , (57) 

L^(t)-L^(0) ) ~t l ' a . an d 



We consider now the stochastic model of dP( a '°\k,t) _ \k\ a p( a ' Hk t) (58) 
Levy flight superdiffusion. dt ' 



VI. FRACTIONAL EQUATION FOR LEVY 
FLIGHT SUPERDIFFUSION 



Taking into account that — \k\ a is the Fourier 

image of the Riesz fractional operator A" ^ 2 = 

d a /d\x\ a , we arrive at the fractional differen- 

T f o ,i t y ,. i ,i tial equation 

if a — 2, the Levy motion becomes the 

Brownian motion with characteristic function „/•„■>/ s i ^ , 

dP(<*>°)(x,t) _ d a P( a >°\x,t) 

p(2-o) (M) = e -^ ? (52 ) Wt ~ d\xf ' ( ' j9) 

obeying the differential equation By means of the direct Fourier transforma- 

tion, one can be convinced of the validity of 
np(2o)/ji f \ two following integral representations of the 

' I ' } = -k 2 p( 2 '°\k,t), (53) Riesz derivative 



d a f(x) = _1_ /■+«> f(x)-f(Q = _L_ ^+oo 2/(a;)-/(a:-0-/(x + 
<9|x| a " K(a)J-oc |x-e| 1+Q " if(a)io 

, (60) 



Here 



Finally, in the case of the asymmetric Levy 
motion, the equation for probability distribu- 
tion becomes 



This equation contains the Feller fractional 
space derivative D^^\ which is determined 
by the relation 



dPM{x,t) _ (aJI) (aJI) 
dt 



x,t). (62) 



Dj*»f(x) 



A(a, j3) 1 + (3 sgn (x - f ) 



/. 



K{a) J-oo \x-£\ 1+a 
MaJ) 2f(x)-(l + (3)f(x- 
K{a) Jo £ 



lf(x) - f(ZM 



l+a 



d£, (63) 



where 

A(a, P) = l + p 2 tg(a%/2). (64) 

A more detailed consideration of fractional 
differential equation for description of Levy 
motion can be found in Refs. [Saichev & Za- 
slavsky, 1997; Uchaikin, 1999, 2000, 2002, 
2003a, 2003b; Uchaikin & Zolotarev, 1999; 
Metzler & Klafter, 2000; Mainardi et ai, 
2001; Metzler & Nonnenmacher, 2002; Lenzi 
et ai, 2003; Gorenflo & Mainardi, 2005; 
Sokolov & Chechkin, 2005; Zaslavsky, 2005]. 

VII. LEVY WHITE NOISE 

Let us come back to the Levy processes. 
The time derivative of the Levy process 



r 



£(t)=dL(t)/dt = L(t) (65) 



is a stationary random process and has anal- 
ogy to the Gaussian white noise, which is the 
time derivative of the Wiener process. The 
Levy process, in fact, is a generalized Wiener 
process. 

Now we derive the characteristic func- 
tional Qt[u] of this Levy delta-correlated 
noise. By definition, we have 



G t [u] = (exp^t J\(r)ar)dr^ = ^exp {i J* u (r) dL (r)| 



^exp |i Km ]T u (0 fc ) [L (r k ) - L (r fc _x)] j j (66) 



= lim ( J] exp {iu (# fc ) [L (r k ) - L (r^)]} ) = lim [T P(u(& k ), At,), 

I 

where ~& k is some internal point of time inter- val (r k -i, r k ) , 5 T = max Ar k , Ar k = r k — Tk-\ 

k 
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(r = 0, r n = t), and P(u('&k), Ar fe ) is the pendence of increments of Levy process L (t). 
characteristic function of increments. To ob- Further from Eqs. (JTJD and ([66]) we obtain 
tain Eq. ([561 we used the statistical inde- 



™ f /•+«>e iu ^ a! -l-zu(7? fc )sina: , 
8* M = Inn || exp < Ar k / - p (x) dx 



k=l 



f n +00 gte^)* _l_ iu (0 ) sina; ] 

exp < hm 2^ Ar fc / = /° ( x ) dx 



k=l 



rt r+oc e tu(r)x — \ — i u ( T \ s i nx ) 

exp \ I dr I - p [x) dx > . 



(67) 



Now we are going to derive a useful func- 
tional formula for the Levy white noise. The 
formula to split the correlation between a 
Gaussian random vector field £ (r, t) and its 
arbitrary functional R[£] was for the first 
time obtained by Furutsu [Furutsu, 1963] and 
Novikov [Novikov, 1965]. For Gaussian ran- 
dom process £ (t) with zero mean it reads 

mm) = / Jf (*' r )(^) dT ' ( 68 ) 

where K(t,r) = (£(X)£(t)) is the correla- 
tion function of Gaussian noise £ (t). We use 



the generalization of Furutsu-Novikov for- 
mula ([651) . obtained by Klyatskin [Klyatskin, 
1974], for arbitrary functional Rt[£] of a non- 
Gaussian random process £ (r), defined on 
the observation interval r £ (0, t), 



iu (t) 



. s 

i5z 



(Rt[t + z\). 

(69) 

Here z (t) is arbitrary deterministic function, 
and $t [u] = In t [iz] . From Eq. ( [671) we ob- 
tain the following expression for variational 
operator in Eq. ([691) 



$t [it] /-+ 00 e iu(i):c - I - iu (t) sin x 



m (t) 



iit (t) x 2 



p (x) dx 



+°° p (x) 



-oo a; jo 



dx / [e iu(t)!/ -cosy]dy. (70) 



r 



Substituting this equation in Eq. ([69]) we ar- rive at 
l 



(t{t)Rt [t + z]) 



+°° p (x) 



dx 



x A 



exp < y 



5z(t) 



— cos y 



{R t [t + z\)dy. 



(71) 



By inserting the operator of functional differ- entiation into the average in Eq. ([71]) and by 
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putting z = 0, we get finally 



p(x) 



dx 



— oo X 



dy. (72) 



VIII. DERIVATION OF KOLMOGOROV'S of the random process X(t) 
EQUATION 



Let us consider now the Langevin equation 
with the Levy white noise source 

X = f(X,t)+g{X,t)Z{t). (73) 



P(x,t) = (6(x-X(t))), (74) 



By differentiating with respect to time the 

well-known expression for probability density and taking into account Eq. (1731 . we obtain 



% = & ^ " ^ 3(x, t) <£ (t) 5(x-X (t))) 



(75) 



r 



By using functional differentiation rules Ref. [Hanggi, 1978], from Eq. (!73|) we get 
and following the same procedure used in 



^S(x-X(t)) = ~g(x,t)S(x-X(t)). 



(76) 



Thus, the variational operator 5/5£ (t) with Eqs. (T72|) . (J75J), and flTBI) . we obtain, after in- 

respect to the function 8 (x — X (t)) is equiv- tegration, the following Kolmogorov's equa- 

alent to the ordinary differential opera- tion for nonlinear system (1731) driven by Levy 

tor —d/dx(g(x,t)). Taking into account white noise [Dubkov & Spagnolo, 2005] 



dP 
~dt 



d[f(x,t)P] + f+°°p(z) 



dx 



exp 



9 , ^' 

-z—g{x,t) 



■ 9 , , 
1 + sin z——q (x, t) 

dx 



dzP. (77) 



We analyze further some different kernel Kolmogorov's equation (ITT]) , related to dif- 
functions p (x) to obtain particular cases of ferent non-Gaussian white noise sources. 
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(a) As a first simple case we consider a 
Gaussian white noise The correspond- 

ing kernel function is p (x) = 2DS (x), where 
D is the noise intensity. After substituting 
this kernel in Eq. (177)) . we obtain the ordi- 
nary Fokker-Planck equation 

9P 9 P D 9 9 P 78 

dt dx dx ^ dx 



(b) For additive driving noise £(t), 
g(X,t) = 1 in Eq. ( 1731) . and the exponen- 
tial operator in Eq. (I77|) reduces to the space 
shift operator. As a result, we find 



dP d r +oc 



d 



P{z) 



P (x — z,t) — P (x, t) + sin z 



dP (x, t) 
dx 



dz. 



(79) 



Equation (1791 is similar to the Kolmogorov- 
Feller equation for purely discontinuous 
Markovian processes [Saichev & Zaslavsky, 
1997; Kamihska & Srokowski, 2004; Dubkov 
& Spagnolo, 2005] 

dP r +OD 

— - = v J w (x — z) P (z,t) dz—vP [x, t) , 

(80) 

where w (x) is the probability density of 
jumps step, and v is the constant mean rate 
of jumps. By putting / (A, t) = 0, omitting 



the term with sinz, and comparing Eq. (|79p 
with Eq. ( IHUl) we find the kernel function for 
such a case 



p (x) = ux 2 w (x) . 



(81) 



For non-Gaussian additive driving force 
^ a - ) (t), with symmetric a-stable Levy dis- 
tribution, the kernel function reads p (x) = 
Q |x| 1_a . As a result, Eq. (1791 takes the fol- 
lowing form 



dP 
~dt 



-I i/(*.«)^+«£ 



P{z,t) -P{x,t) 



\x 



\a+l 



dz 



(82) 



r 



and describes the anomalous diffusion in form 
of symmetric Levy flights. 

In accordance with the definition of Riesz 
derivative (60), Eq. (82) can be written as 

dP 
~dt 



d d a P 
-^[f^t)P] + D^— 
ox a \x\ 



3) 



where (see Eq. (61)) 
D = K (a)Q 



r(a + 1) sin (vra/2) ' ^ 

For the first time, the fractional Fokker- 
Planck equation (|83j) for Levy flights in the 



potential profile U (x) (with —U' = f(x,t) 
was obtained directly from Langevin equa- 
tion 



X = -U' (A) + £ (Q) (t) 



(85) 



by replacing the white Gaussian noise £ (t) = 
£,( 2 ' (t) with the symmetric Levy a-stable 
noise £ (a) (t), in Refs. [Ditlevsen, 1999b; 
Yanovsky et ai, 2000; Garbaczewski & 
Olkiewicz, 2000] (see also [Schertzer et ai, 
2001]). However, some attempts were under- 
taken before in Refs. [Fogedby, 1994a, 1994b, 
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1998; Jespersen et ai, 1999]. Recently us- 
ing a different approach it was derived in 
[Dubkov & Spagnolo, 2005]. 

IX. STATIONARY PROBABILITY 
DISTRIBUTIONS FOR LEVY FLIGHTS 



Let us consider, as in Ref. [Chechkin et ai, 
2002a], the symmetric smooth monostable 
potential U (x) = ^x 2m / (2m) (m = 1,2,...). 
The Eq. (|89|) . therefore, transforms into the 
following differential equation of (2m — 1)- 
order 



First of all, we can try to evaluate the 
stationary probability distribution P st (x) of 
Levy flights in the potential profile U(x) from 
Eq. (|83|) . Of course, this evaluation is impos- 
sible for any potential profile, but the poten- 
tial U(x) should satisfy some constraints. It 
is better to apply Fourier transform to the 
integro-differential equation (J83l) and to write 
the equation for the characteristic function 



P(k,t) 



+oo 



e ikx P (x, t) dx. 



(86) 

After simple manipulations we find (see 
Eq(58)) 



dP 
~dt 



—— = —ik 



P. 



(2m-l) 



si 



■l) m+1 /? 2m_1 



ia-1 



sgn k-P st = , 

(90) 

where (3 = 2m ~^j 1 D Jj. As it was proved by 
analysis of Eq. ( 1901) for small arguments k 



in Ref. [Chechkin et al., 2002a], the station- 
ary probability distribution P st (x) has non- 
unimodal shape and power tails 



Pst (X) 



1 



i2m+a-l ' 



X 



oo . 



\x\ 



(91) 



According to Eq. (I9TI) , we have a confinement 
of Levy flights (finite variance of particle's 
coordinate) in the case when 



2m > 4 



a. 



(92) 



e tkx U'(x)P (x, t) dx-D \k\ a P. 

(87) 

For smooth potential profiles U (x), expand- 
ing in power series in a neighborhood of the 
point x = 0, we can rewrite Eq. (87) in the 
operator form 



OP 
~dt 



-ikU' 



'dk 



P-D\k\ a P. 



Because of the Levy index a G (0,2], a 
confinement takes place for all values of a, 
starting from quartic potential (m = 2). Ex- 
act solution of Eq. (1911 can be only obtained 
for the case of Cauchy noise source £W (t) 
(a = 1). Due to the symmetry of the char- 
acteristic function P st (—k) = P st (k), we can 
reduce Eq. f )90l) to a linear differential equa- 
tion with constant parameters 



In particular, for stationary characteristic 
function, from Eq. flHBT) we get 




zD|fc| a_1 sgnfc-P st = 



39) 



,(2m-l) 



st 



(fc>0). 

(93) 

From the corresponding characteristic equa- 
tion 



A 



2m- 1 



m n2m-l 



(94) 



where sgn/c is the sign function. Unfortu- we select the roots with negative real part, 
nately, one cannot solve analytically Eq. (|89|) which are meaningful from physical point of 
for arbitrary potential U (x) and arbitrary view. The general solution of Eq. ([93]), there- 
Levy exponent a. fore, reads 
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[(m " 1)/2] f 7T (m - 21 - 1) 

Pst(k)= Aiexp < -(3 \k\ cos' 

i=o I 



2m 



( . 7r(m-2Z-l) 
cos I p I ft I sm ^ - ipi 



(95) 



where the quadratic brackets in the upper where j = 1,2, . . . ,m — 1. Now substituting 
limit of the sum [(m — l)/2] denote the in- Eq. (195]) in Eq. (l96l) we obtain 
teger part of the enclosed expression. The 
unknown constants Ai and tpi can be calcu- 
lated from the conditions 



Pst (0) 



P. 



(2j-l) 



+0) = o, 



(96) 



[(m-l)/2] 



«=0 



A/ cos v?/ 



[(m-l)/2] 
1, Ai COS 

Z=0 



7T (2j - 1) (m + 2/) 
2m- 1 



(j = l,2,...,m-l). 

(97) 



Making the reverse Fourier transform in distribution (SPD) of the particle coordinate 
Eq. (1951) we find the stationary probability 



Pst (X) 



[(m-l)/2] 

2=0 



x 2 cos 



7r(m-2Z-l) , 



+ /3 2 COS 



7r(m-2/-l) 
2m- 1 



x 4 - 2x 2 /3 2 cos 2g±ii + /3 4 
I 



(98) 



The parabolic potential profile U (x) = 7x 2 /2 For quartic potential (m = 2), from the set 

corresponds to a linear system (1751) . In this of equations (19~7I) . we find A = 2/ a/3 and 

situation, from Eqs. (19T1) and (I9"g|) we easily <^ = 7r/6. Substituting these parameters in 

obtain the following obvious result Eq. (198]) we obtain 



P,t (x) 



P 



7T (x 2 + (3 2 ) 



(99) 



Pst (X) 



P 3 



(100) 



71 (X 4 - X 2 P 2 + /? 4 ) ' 

that is due to the stability of the Cauchy dis- 
tribution (I99j) . the probabilistic characteris- which coincides, for (3=1, with the result 
tics of driving noise increments (see Eq. (H2~]) ) obtained in Refs. [Chechkin et ai, 2002a, 
and Markovian process X (t) are similar (see 2006]. The plots of stationary probability 
also Ref. [West & Seshadri, 1982]). distributions (11001) for Levy flights in sym- 
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metric quartic potential, for different values 
of the parameter j3, are shown in Fig. 1. 




FIG. 1 Stationary probability distributions for 
Levy flights in symmetric quartic potential 
U(x) = 7x 4 /4 for different values of dimen- 
sionless parameter j3: (1) = 0.5, (2) = 1, 
(3) = 1.5. 



The superdiffusion in the form of Levy 
flight gives rise to a bimodal stationary prob- 
ability distribution, when the particle moves 
in a monostable potential. This bimodal dis- 
tribution is a peculiarity of Levy flights. In 
fact the ordinary diffusion of the Brownian 
motion is characterized by unimodal SPD. 
The SPD of superdiffusion has two maxima 
at the points x = ±/3/\/2, with the value 
{P s t) max = 4/ (3ti0). Since the value of the 
minimum is P st (0) = 1/ (vr/3), the ratio be- 
tween maximum and minimum value is con- 
stant and equal to 4/3. The variance of the 
particle coordinate, obtained from Eq. (TjQ) is 
finite: (X 2 ) st = /3 2 . As a result, the prob- 
ability distribution becomes more wide with 
increasing parameter j3 - 



'D/j, that is with 
decreasing the steepness 7 of the quartic po- 
tential profile, or with increasing the noise 
intensity D. 

A detailed analysis of the solution of the 
differential equation (1901 . for arbitrary Levy 



index a and quartic potential (m = 2) was 
performed in Refs. [Chechkin et ai, 2002a, 
2004]. In Fig. 2 the profiles of SPD (ob- 
tained by an inverse Fourier transformation) 
in symmetric quartic potential are shown for 
the different Levy indices from a = 1, at the 




FIG. 2 Forms of stationary probability distribu- 
tion in the symmetric quartic potential for dif- 
ferent Levy indices, from a = 1 till a = 2. From 
Ref. [Chechkin et ai, 2002a]. 



top of the figure, up to a = 2 at the bot- 
tom [Chechkin et ai, 2002a]. It is seen that 
the bimodality is most strongly expressed for 
a = 1 (Cauchy stable noise source). By in- 
creasing the Levy index, the bimodal profile 
smoothes out, and, finally, it turns into a uni- 
modal one at a = 2, recovering the Boltz- 
mann distribution. 

Carrying out analogous procedure we ob- 
tain the stationary probability distributions 
for the cases m = 3, 4, 5 [Dubkov & Spagnolo, 
2007] 



Pst (x) 



t 

vr {x 2 + (3 2 ) {x 4 - 2fJ 2 x 2 cos vr/5 + /3 4 ) ' 
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Pst (X) 
Pst (X) 



£ 

7i (x A - 2p 2 x 2 cos tt/7 + p A ) (x 4 + 2(3 2 x 2 cos 2tt/7 + p 4 ) ' 

£ 

n (x 2 + p 2 ) (x 4 - 2p 2 x 2 cos vr/9 + p 4 ) (x 4 + 2/3 2 x 2 cos 4vr/9 + /3 4 



(101) 



The plots of distributions (11011) . for dif- 
ferent values of parameter /3, are respectively 
shown in Figs. 3-5. 

It must be emphasized that according to 
Figs. 3-5, these distributions remain bimodal 
and have the same tendency with increas- 
ing j3, but the ratio between maximum and 
minimum increases with increasing m. From 
Eqs. (HOOP and (jlOip we see that the second 
moment of the particle coordinate is finite for 
m > 2, which confirms the inequality ( |92l) . 
This means that there is a confinement of 
the particle motion due to the steep potential 
profile, even if the particle moves according 
to a super diffusion in the form of Levy flights. 
The presence of two maxima is a peculiarity 
of the super diffusion motion. Because of the 
fast diffusion due to Levy flights, the parti- 
cle reaches very quickly regions near the po- 
tential walls on the left or on the right with 
respect to the origin x = 0. Then the par- 
ticle diffuses around this position, until a 
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FIG. 3 Stationary probability distributions for 
Levy flights in symmetric potential U(x) = 
7j; 6 /6 for different values of dimensionless pa- 
rameter (3: (1) [3 = 0.5, (2) & = 1, (3) = 1.5. 



new flight moves it in the opposite direction 




-2-1012 



FIG. 4 Stationary probability distributions for 
Levy flights in symmetric potential U(x) = 
7x 8 /8 for different values of dimensionless pa- 
rameter (3: (1) p = 0.5, (2) p = I, (3) p = 1.5. 

to reach the other potential wall. As a re- 
sult, the particle spends a large time in some 
symmetric areas with respect to the point 
x = 0, differently from the Brownian diffu- 
sion in monostable potential profiles. These 
symmetric areas lie near the maxima of the 
bimodal SPD. For fixed D and m, these max- 
ima are closer or far away the point x = de- 
pending on the greater or smaller steepness 7 
of the potential profile. This corresponds to a 
greater or smaller confinement of the particle 
motion. Of course, such confinement is more 
pronounced for greater m, that is for steeper 
potential profiles. 

On the basis of Eqs. (|99|)-(fTQll and the 
known behavior of density tails (1911) . we can 
write the general expressions for stationary 
probability distribution in the case of po- 
tential U (x) = 7X 2 " 1 / (2m) with odd m = 
2n + l [Dubkov & Spagnolo, 2007] 
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FIG. 5 Stationary probability distributions for 
Levy flights in symmetric potential U(x) = 
7j; 10 /10 for different values of dimensionless pa- 
rameter (3: (1) = 0.5, (2) = 1, (3) p = 1.5. 



^ ^ = 7T (X 2 + /5 2 ) 2 X 4 _ 2/ g2 x 2 CQS ^ ( 4/ + J) / ( 4n + !)] + ^4 > ( 102 ) 

I 



and even m = 2n 



g4n—l n—1 

^ ^ = "V" 2 X A _ 2fj2 x 2 cos [tt (4/ + 1) / ( 4n _ 1)] + ^4 • ( 103 ) 
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The strong proof of non-unimodality of Because of the symmetry of the SPD P st (x) , 

the SPD for symmetric monostable potential Eq. (11051) gives 
U (x) = \x\ c /c in the case c > 2 was given 

in Ref. [Chechkin et al, 2004]. Indeed, from roo p ( z \ _ p (q\ 

Eq. (JE2D we have for SPD J - 1 '— x - 1 ' = 0. (106) 



d 

dx 



I c — 1 n 

x\ sgn x ■ K 



St 



d a P t F° r unimodal probability distribution with 



^ d\x\ a ~ ^ the maximum at the origin, the integral in 

the left side of Eq. (1 106f) should be negative, 

As a result, from Eqs. pD and (60) at the which contradicts Eq . pg| , 

point x = we obtain The estimation of bifurcation time for 

transition from unimodal initial distribution 

+ °° -Pst ~ fit (0) _ q to bimodal stationary one for the quartic po- 

-oo tential (c = 4) was done in Refs. [Chechkin 

(105) et al, 2004, 2006]. The dependence of this 



d a P< 



si 



d \x\ 



2 = 
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bifurcation time t\2 from Levy index a is 
plotted in Fig. 6. Also authors proved an ex- 




a. s 



o ■ 

l.o 1.5 a.o 

a 

FIG. 6 Bifurcation time t±2 versus Levy index 
a for quartic potential. Black dots: bifurca- 
tion time deduced from the numerical solution of 
the fractional Fokker-Planck equation. Dashed 
and solid lines: two subsequent approximations. 
From Refs. [Chechkin et al, 2004, 2006]. 

istence of a transient trimodal state between 
initial unimodal and final bimodal ones. This 
evolution, shown in Fig. 7, can be only ob- 
served for monostable potential with c > 4 
and for fixed values of the Levy index a. 
The corresponding bifurcation times of tran- 
sitions ti3 (unimodal — * trimodal) and t 32 
(trimodal — ■> bimodal) versus Levy index a, 
with potential exponent c = 5.5, are plotted 
in Fig. 8. 

X. BARRIER CROSSING 

The problem of escape from metastable 
states investigated by Kramers [Kramers, 
1940] is ubiquitous in almost all scientific 
areas [Hanggi et al, 1990; Spagnolo et al., 
2007]. Since many stochastic processes do 
not obey the Central Limit Theorem, the cor- 
responding Kramers escape behavior will dif- 
fer. An interesting example is given by the te- 
stable noise-induced barrier crossing in long 
paleoclimatic time series [Ditlevsen, 1999a]. 
Another new application is the escape from 
traps in optical or plasma systems [Fajans & 




FIG. 7 The evolution of the probability dis- 
tribution for a = 1.2 and c = 5.5 from 
unimodal through trimodal to bimodal. From 
Refs. [Chechkin et al, 2004, 2006]. 
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FIG. 8 Bifurcation times t±s and £32 versus Levy 
index a for the power potential with the exponent 
c = 5.5. From Ref. [Chechkin et al., 2004]. 

Schmidt, 2004]. 

The main tools to investigate the barrier 
crossing problem for Levy flights are the first 
passage times, crossing times, arrival time 
and residence times. We should emphasize 
that the problem of mean first passage time 
(MFPT) meets with some difficulties because 
free Levy flights represent a special class of 
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discontinuous Markovian processes with infi- 
nite mean squared displacement. First of all, 
the fractional Fokker-Planck equation ( |83l) is 
integro-differential, and the conditions at ab- 
sorbing and reflecting boundaries differ from 
the usual conditions for ordinary diffusion. 
Superdiffusion motion is characterized by the 
presence of jumps, and, result, a parti- 
cle can reach instantaneously the boundary 
from arbitrary position. One can mention 
some erroneous results for Levy flights ob- 
tained in Ref . [Gitterman, 2000] , because au- 
thor used the traditional conditions at two 
absorbing boundaries (see the related corre- 
spondence [Yuste & Lindenberg, 2004; Git- 
terman, 2004]). The numerical results for the 
first passage time of free Levy flights con- 
fined in a finite interval were presented in 
Ref. [Dybiec et al, 2006]. The complexity 
of the first passage time statistics (mean first 
passage time, cumulative first passage time 
distribution) was elucidated together with a 
discussion of the proper setup of correspond- 
ing boundary conditions, that correctly yield 
the statistics of first passages for these non- 
Gaussian noises. In particular, it has been 
demonstrated by numerical studies that the 
use of the local boundary condition of vanish- 
ing probability flux in the case of reflection, 
and vanishing probability in the case of ab- 
sorbtion, valid for normal Brownian motion, 
no longer apply for Levy flights. This in turn 
requires the use of nonlocal boundary condi- 
tions. Dybiec with co-authors found a non- 
monotonic behavior of the MFPT for two ab- 
sorbing boundaries, with the maximum being 
assumed for a = 1 (see Fig. 9), in contrast 
with a monotonic increase for reflecting and 
absorbing boundaries. 

According to the Kramers law, the proba- 
bility distribution of the escape time from a 
potential well with the barrier of height Uq, 
has the exponential form 



p{t) 



1 



exp 



t 



(107) 




FIG. 9 Mean first passage time versus Levy in- 
dex a of confined motion between two absorb- 
ing boundaries driven by stable symmetric Levy 
white noise. From Ref. [Dybiec et al., 2006]. 



T c = C exp 



Uo 
D 



(10£ 



where C is some positive prefactor and D 
is the noise intensity. The problem how 
the stable nature of Levy flight processes 
generalizes the barrier crossing behavior of 
the classical Kramers problem was investi- 
gated, both numerically and analytically, in 
Ref. [Chechkin et al, 2003b, 2005, 2006, 
2007]. Authors considered Levy flights in a 
bistable potential U (x) by numerical solu- 
tion of Eq. (1851) . It was shown that although 
the survival probability decays again expo- 
nentially as in Eq. (11071) . the mean escape 
time T c has a power-law dependence on the 
noise intensity D 



C(a) 



(109) 



with mean crossing time 



where the prefactor C and the exponent \x de- 
pend on the Levy index a. Using the Fourier 
transform, i.e. Eq. (IHHl) . the mean escape 
rate was found for large values of 1/D in 
the case of Cauchy stable noise (a = 1) in 



22 



the framework of the constant flux approx- 
imation across the barrier. The probability 
law and the mean value of escape time from 
a potential well for all values of the stabil- 
ity index a G (0,2), in the limit of small 
Levy driving noise, were also determined in 
the paper [Imkeller & Pavlyukevich, 2006] by 
purely probabilistic methods. Escape times 
had the same exponential distribution (11071) . 
and the mean value depends on the noise in- 
tensity D in accordance with Eq. (11091) with 
/i(a) = 1 and pre-factor C depending on a 
and the distance between the local extrema 
of the potential. 

The barrier crossing of a particle driven 
by symmetric Levy noise of index a and in- 
tensity D for three different generic types 
of potentials was numerically investigated in 
Ref. [Chechkin et al, 2007]. Specifically: (i) 
a bistable potential, (ii) a metastable poten- 
tial, and (iii) a truncated harmonic potential, 
were considered. For the low noise intensity 
regime, the result of Eq. (11091) was recovered. 
As it was shown, the exponent n(a) remains 
approximately constant, fi ~ 1 for < a < 2; 
at a = 2 the power-law form of T c changes 
into the exponential dependence (I108p . It 
exhibits a divergence-like behavior as a ap- 
proaches 2. In this regime a monotonous in- 
crease of the escape time T c with increasing a 
(keeping the noise intensity D constant) was 
observed. For low noise intensities the escape 
times correspond to barrier crossing by multi- 
ple Levy steps. For high noise intensities, the 
escape time curves collapse for all values of a. 
At intermediate noise intensities, the escape 
time exhibits non-monotonic dependence on 
the index a as in Fig. 9, while still retains the 
exponential form of the escape time density. 

The first arrival time is an appropriate pa- 
rameter to analyze the barrier crossing prob- 
lem for Levy flights. If we insert in frac- 
tional Fokker-Planck equation (183|) a 5-sink of 
strength q (t) in the origin, we obtain the fol- 
lowing equation for the non-normalized prob- 
ability density function P (x, t) 



dP d <9 a P 

_ = _ [W]+D ^- 5 ( t )* w , 

(110) 

from which by integration over all space we 
may define the quantity 



?(*) 



d 
dt 



P (x, t) dx, 



111) 



which is the negative time derivative of the 
survival probability. According to definition 
(II 111) , q (t) represents the probability den- 
sity function of the first arrival time: once 
a random walker arrives at the sink it is 
annihilated. As it was shown in the paper 
[Chechkin et al., 2003b] for free Levy flights 
(U (x) = 0), the first arrival time distribution 
has a heavy tail 



q (t) ~ t 1/a ~ 



(112) 



with exponent depending on Levy index 
a (1 < a < 2) and differing from universal 
Sparre Andersen result [Sparre Andersen, 
1953, 1954] for the probability density func- 
tion of first passage time for arbitrary Marko- 
vian process 



p(t)~f 



-3/2 



(113) 



In the Gaussian case (a = 2), the quantity 
(11121) is equivalent to the first passage time 
probability density (I113p . From a random 
walk perspective, this is due to the fact that 
individual steps are of the same increment, 
and the jump length statistics therefore en- 
sures that the walker cannot hop across the 
sink in a long jump without actually hitting 
the sink and being absorbed. This behavior 
becomes drastically different for Levy jump 
length statistics: there, the particle can eas- 
ily cross the sink in a long jump. Thus, before 
eventually being absorbed, it can pass by the 
sink location numerous times, and therefore 
the statistics of the first arrival will be differ- 
ent from that of the first passage. The result 
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(11131) for Levy flights was also confirmed nu- 
merically in the paper [Koren et at, 2007]. 

At last, the nonlinear relaxation time tech- 
nique is also suitable for investigations of 
Levy flights temporal characteristics. Ac- 
cording to definition, the mean residence time 
in the interval (Li, L 2 ) reads 

T(x ) = dt P(x,t\x Q ,0)dx, (114) 

JO JLi 

where Xq is the initial position of all parti- 
cles (x G (Li,L 2 )) and P(x, t\x o ,0) is the 
probability density of transitions. We do not 
need to think about the boundary conditions 
in this case because we are concerned with 
the overall time spent by the particle in the 
fixed interval. Changing the order of integra- 
tion in Eq. (11141) we obtain 

T (x ) = [ L2 Y (x, x , 0) dx, (115) 

where Y (x, xo, s) is the Laplace trans- 
form of the transient probability density 
P (x, t\ xo, 0) 

YCx,x ,s)— / P(x,t\ x ,0)e~ st dt. 
Jo 

(116) 

Making the Laplace transform in the frac- 
tional Fokker-Planck equation fl83l) and 
taking into account the initial condition 

P (x, 0| xo, 0) = 5 (x — xo), we get 

^ [V (x) Y] + D-^s -sY = -5{x- x ) . 

(117) 

If we put s = in Eq. fl 1 1 T|) and make the 
Fourier transform we obtain 

/ d \ ~ 1 ~ e ikx ° 

U' l-i — j Y-iD {kr 1 sgn (k) Y=—, 

(118) 

where 

Y(k,x )= Y (x, x , 0) e +tkx dx. (119) 

J — oo 



After solving Eq. (11181) we can calculate the 
mean residence time as (see Eqs. (11151) and 

(HID) 

l r+oc e —ikLi _ e ~ikL2 _ 

T ( x o) = 7T- ; Y(k,x )dk. 

LlXl J-oo K 

(120) 

Equations fll 18j) and (I120p are useful tools to 
analyze the temporal characteristics of Levy 
flights in different potential profiles U(x). 

XI. CONCLUSIONS 

In this tutorial paper, after some short 
historical notes on normal diffusion and su- 
perdiffusion, we introduce the Levy flights as 
self-similar Levy processes. After the def- 
inition of the strictly stable random vari- 
ables, the subfamily of the Levy motion 
is introduced with the fractional differen- 
tial equation for Levy flight superdiffusion. 
We used then functional analysis approach 
to derive the fractional Fokker-Planck equa- 
tion directly from Langevin equation with 
symmetric a-stable Levy noise. This ap- 
proach allows to describe anomalous diffu- 
sion in the form of Levy flights. We obtained 
the general formula for stationary probabil- 
ity distribution of superdiffusion in symmet- 
ric smooth monostable potential for Cauchy 
driving noise. All distributions have bimodal 
shape and become more narrow with increas- 
ing steepness of the potential or with de- 
creasing noise intensity. We found that the 
variance of the particle coordinate is finite 
for quartic potential profile and for steeper 
potential profiles, that is a confinement of 
the particle in a superdiffusion motion in the 
form of Levy flights. As a result, we can eval- 
uate the power spectral density of a station- 
ary motion. We have also discussed recently 
obtained analytical and numerical results for 
time characteristics of Levy flights. Special 
attention was given for some difficulties with 
formulation of the correct boundary condi- 
tions for mean first passage time problem. 
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As it was shown, the arrival and residence 
times are more appropriate characteristics for 
investigations of Levy flights in different po- 
tential profiles. 
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